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1.0 


INTRODUCTION 


The NASA Lewis Research Center employs a general computer program NNEP (Ref.l) 
for calculating the thermodynamic performance of jet propulsion engines. To calculate 
off-design engine performance, the NNEP user must input component maps defining the 
characteristics of the various components over their full range of operating con- 
ditions . 

For early cycle analysis of advanced propulsion systems, these map characteris- 
tics are not generally known because the geometry of the component has not been 
specified. Furthermore, the typical user of NNEP is not sufficiently knowledge- 
able and/or cannot afford the time to do an extensive design followed by an off-design 
analysis of the component in question to define the map characteristics. Typically, 
in this early stage, the user scales some available map. 

The available methods for scaling maps can lead to significant errors in com- 
ponent representations. A traditional method of scaling a compressor map retains the 
flow speed relation of the base map and applies a constant pressure rise scalar cal- 
culated at the design point. Size scaling is, of course, legitimate. The accuracy of 
such a procedure can be considerably improved by using parametrically generated com- 
ponent maps. A parametric component representation can be thought of as a scaling pro- 
cedure which uses the key design point parameters to impact the fundamental differences 
in the map characteristics when generating the component maps. 

The objective of the present study is an improved method of representing single 
stage flow modulating fans and/or compressors of either the axial or centrifugal types 
when performing calculations of off-design performance for advanced air-breathing jet 
engines. The axial flow machines include flow modulation through either variable inlet 
guide vanes (VI GV) or variable rotor pitch (VPF) . The centrifugal compressors include 
only the variable inlet guide vane (VIGV) option. This method, which is a computer pro- 
gram called MODFAN (i.e. , flow modulating fan) is compatible in both form and format 
with the cycle program of Reference 1 and the example map representation of Reference 
2 . 


This report is a user’s manual for the computer program MODFAN and contains a 
description of the input-output data, values of typical inputs, and sample cases. A 
brief description of the engineering analysis used to generate a program is given near 
the end of the report. 1 



The program uses design point data and semi-empirical correlations to 
generate off-design values of corrected f low, efficiency , and pressure ratio over a 
range of blade angles, corrected speeds, and pressure ratio parameters specified 
by the user. 
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2 . 0 PROGRAM STRUCTURE 


A flow chart showing the flow of control in the MODFAN program is shown in 
Figure 1. The program first displays a description of the design point variables 
together with the default values. The user can then change the default values as 
desired. The input is then checked and the updated input displayed. This completes 
the input phase of the program. The design point calculations are then carried out. 
These calculations determine the blade row overall geometries from the design 
point input. 

Once the design point calculations have been completed, the off-design calcula- 
tion can begin. The calculations for each of the desired angles, speeds, and *R values 
are carried out in a set of nested DO loops. For each angle and speed, the min-loss 
point (i.e., the map backbone) is first located. Then the values of work coefficient 
both at stall and at a pressure ratio of unity are determined. These values of work 

coefficient form the upper and lower limits of the speed line. The R values are then 

converted into equivalent work coefficient values, and the fan characteristics 
determined. When the inner or R DO loop is completed the speed is incremented and the 

calculation repeated: when all speeds are completed, the angle is incremented. 

After the off-design calculation has been completed, MODFAN writes the three 
data sets required for the computer program NNEP input. 


*R is a measure of the relative distance along a speed line starting from R=1 at 
stall. 
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READ INPUT 
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FIGURE 1. FLOW CHART SHOWING FLOW OF CONTROL IN MODFAN 






3.0 PROGRAM INPUTS 


Most of the MODFAN inputs are of the free-field format (NAMELIST) type, and 
begin in column two. The only exception is the initial type switch which requires 
either a 1 or 0 as a response. The program first gives a brief description of the 
variables used in the NAMELIST INPUT. The default settings of these variables are 
then displayed. The user can then enter changes in the design point values and/or 
the range of angle settings, speeds, and R values desired. If the user wishes the 
program to generate a value for the design point efficiency, a zero value should be 
input for EFFD. The program will then echo the updated NAMELIST and go into exe- 
cution. Upon completion, the program will display a message to the effect that the 
NASA output files have been written on file codes 30, 31, and 32. 

A sample of this showing the terminal conversation is shown in Figure 2 . 

This example is for a axial flow fan having fixed geometry. Note that all of the 
design point values have been reset as well as the range of inlet guide vane angles 
and corrected speeds. The range of R values has not been altered. The R values are 
used to fix a point on a speed line. The R value is unity at the stall line and in- 
creases along a constant speed line as the flow increases. The algorithm used in 
MODFAN forces a value of R equal to two at the min-loss point which is slightly 
below the peak efficiency on the speed line. The concept of min-loss is discussed 
in detail in Reference 3 . 

The fan map generated by the program from the input values given in Figure 2 
is shown in Figure 3. The locus of the R=1 and R=2 lines have been indicated on 
the figure. 
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0R10W U !\V3E IS 

OF POOR QUA li ! v 


PARAMETRIC FLOW MODULATING FAN 


*M*0*D*F*A*N“ 


AXIAL OR CENT COMPRESSOR? 

-0 


1 *CENT t0 OR CR-AXIAL 


DESCRIPTION OF INPUT VARIABLES IN MODFAN PROGRAM 
(NAMELIST INPUT) 


ITYPE-1 VIGV 

ITYPE-2 VPF ( NO IGV ) 

DESIGN POINT VALUES OF; 


PQPD PRESSURE RATIO! TOTAL -TO-TOTAL ) 

EFFD AO IABATIC EF F IC I ENC Y( TOTAL-TO TOTAL) 

W1D INLET CORRECTED FLOW 

PSID ROTOR EXIT MERIDIONAL LOADING! 2G0*DH/U2**2 ) 

WQA1D INLET CORRECTED FLOW PER UNIT AREA 

XMZ2D MERIDIONAL MACH NO AT ROTOR EXIT 

SMO CONSTANT SPEED STALL MARGIN! PERCENT ) 

STP1D IGV METAL ANGLE 

GEOMETRY SPECIFICATIONS 


A3QA2 STATOR INLET TO ROTOR EXIT AREA RATIO 

RRAT INLET HUB TO TIP RADIUS RATIO 

ROOI RATIO OF ROTOR EXIT TO INLET RADII (MERIDIONAL LINE) 

MAP RANGE (ARRAY VALUES MUST BE IN ASCENDING ORDER) 


NR 


NO OF R 

VALUES TO BE CALCULATED 

AR 


ARRAY OF 

R VALUES! I “STALL ,2“MIN-LOSS 

NSPDS 


NO OF SPEED LINES 

APCNC 


ARRAY OF 

SPEEDS (DECIMAL PERCENT) 

NSTP1 


NO OF IGV SETTINGS 

ASTP1 


ARRAY OF 

IGV ANGLES! -10 TO 40 DEG) 

NROT1 


NO OF ROTOR PITCH SETTINGS 

AROT1 


ARRAY OF 

ROTOR PITCH ANGLES! -5 TO 5 

NAMELIST 

INPUT 


I TYPE 

a 

I, 


PQPD 

a 

l .6550, 

EFFD - 0.0340, 

WID 

a 

223.8000, 

PSID ■ 0.0226, 

VQAiD 

a 

41 .0700, 

XM22D - 0.4S00, 

SMD 

a 

14.1 000 , 

STP1D - 0. 

A3QA2 

a 

0.9600, 

RRAT - 0.5000, 

ROQI 

a 

1 . 0000 , 


NR 

a 

1 1 


AR( I ) 


a 


l 


1 .0000, 

1.2000, 1.4000, 

5 


1 .8000, 

2.0000. 2.2000, 

9 


2 . 6000 , 

2.0000, 3.0000, 


1 . 6000 , 
2 . 4000 , 


Figure 2. MODFAN Sample Terminal Conversation. 
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ORIGINAL PAtfi £ 
OF POOR 0<>^ t J 


NSPDS - 6 

APCNC< I ) - 

1 0.5 000, 0 . 7001 0, 0 . 8000 , 

5 1.0000, 1.1000, 

NSTP1 - 3 

ASTPKI) - 

1 0 . , 20.0000, 40.0000, 

NROT1 - 0 

AROTK I > - 

1 0 . 

END NAMELIST INPUT 

ENTER CHANGES TO NAMELIST INPUT 
-SINPUT NSTP1/ASTP1 *0.0, 

-PQPD-1 .62,EFFD-.84 ,VQA1D-40.0,VIO-219.2 ,PSIO- . 863, 


0.9000, 


-XMZ2D".5,SMD-20.0,STP1D«0.0,A3QA2-1 .05 , RRAT- . 5 , ROQI- 1 .0, 
-NSPDS/APCNC-.5, .7, .9, 1 .0, 1 . 1 , 


-$ 


NAMELIST INPUT 

ITYPE - 1, 

POPD - 1.6200, EFFD - 0.8400, 

V1D - 219.2000, PSID - 0.8630, 

WGA1D - 40.0000, XMZ2D - 0.5000, 

SMD - 20.0000, STP10 » 0. 

A3QA2 - 1.0500, RRAT - 0.5000, 

ROQI - 1.0000, 

NR - 11 

AR( I ) 

1 1.0000, 1.2000, 1.4000, 1.6000, 

5 1.8000, 2.0000, 2.2000, 2.4000, 

9 2.6000, 2.8000, 3.0000, 

NSPDS - 5 

APCNC ( I > * 

1 0.5000, 0.7000, 0.9000, 1.0000, 

5 1.1000, 

NSTP1 * I 

ASTP 1 ( I ) * 

1 0 . 

NROT1 - 0 

AROTK I) » 

1 0 . 

END NAMELIST INPUT 

OUTPUT ON IFC-30,31,32 


PROGRAM COMPLETE & 


Figure 2. MODFAN Sample Terminal Conversation (Continued). 
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4 . 0 PROGRAM OUTPUTS 


The basic output from the program consists of three tables. These tables 
show the variation in corrected flow, efficiency, and pressure ratio for each of 
the vane or rotor angles, R-values, and corrected speeds specified in the input. 

The output tables for the default cases are shown on pages 15 thru 17 and 20 thru 
22. The table structure is compatible with NASA cycle deck requirements given in 
Reference 2 (pages 23 and 24) . 

The output tables can be visualized as three dimensional, composed of a 
series of planes with each plane assigned a value of angle position, STP1 or R0T1. 
Then in each angle plane, the dependent variable (ordinate axis) is a function of 
corrected speed, SPED, and R value. The dependent variables are respectively 
corrected flow, total-to-total efficiency, and pressure ratio. 

For example, in the output table on page 20 the 12 lines of the dependent 
variable correspond to the 6 values of corrected speed, two lines per speed. Within 
each speed there are 11 R values. 

It should be noted that for pressure ratios less than unity the efficiency is 
negative. These efficiency values are not incorrect, however, the efficiency be- 
havior in this region makes curve fitting and interpolation of efficiency values 
extremely difficult. For this reason many engine manufactures use some form of locus 
or temperature ratio parameter rather than efficiency for interpolation. The solution 
used here was to simply discard the information below unity pressure ratio and to 
display the solution for unity pressure ratio for all values of R at which the 
pressure ratio is less than unity. This means that identical values of pressure ratio 
(PR=1.0), flow and efficiency (EFF=0.0) will appear in the output table on any speed 
line where the value of R results in a pressure ratio less than unity. 
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5.0 PROGRAM DIAGNOSTICS 


The MODFAN computer program contains error printouts to aid the user in 
trouble shooting. A list of the error messages and their meaning are given 
below. 


ERROR IN SUBROUTINE MINLOS. 

PCN - F7.3, STP1 « F7.3, ROTl - F7.3, ERR = F7.4 

Failure to find the min-loss or "backbone location" 
at the input corrected speed and angle. The rotor 
is not choked. 

ERROR IN SUBROUTINE ROTORC 

TANA1 = F7.3, PCN = F7.3, PSI2J = F7.3, ERR = F7.4 

Failure in the Newton-Raphson loop which solve the 
continuity equation across the rotor for the case 
where the rotor is choked. 

NO STALL INTERSECTION FOUND: PCN = F7.2 

Failure in the STALLX Subroutine. No stall intersection 
found at the specified speed and *J>>2. 
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There are three iterations in the program which are balanced using the Method 
of False Position. This method is contained in the subroutine QIREXX. A maximum of 
25 passes is allowed for any single iteration to balance. If the iteration d 0 »s not 
balance within the specified tolerance the error message shown below will printout 
with the number of the offending iteration in the 15 format field. 

QIRE CTR ERROR (CALLING LINE = 15) 


QIRE 

CALLING 


LOOP 

ROUTINE 

COMMENTS 

1 

MINLOS 

Calculation of Min-Loss or "Backbone" 
location at the input corrected speed 
and angle with Rotor choked. 

2 

STALLX 

Calculation of the intersection of 
speed line and stall line (R=l) . 

3 

PRCHOP 

Calculation of the corrected flow at 
unity pressure ratio. 
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6 . 0 EXAMPLE CASES 


Two example cases are given in order to illustrate the use of the program. 
The first case utilizes the default settings to generate a single stage fan having 
variable inlet guide vanes (VIGV) , The second case uses the default settings to 
generate a single stage fixed geometry centrifugal compressor. 

A complete record of the two terminal sessions including a listing of the 
output tables is given on the following pages. The program inputs and outputs have 
been discussed previously in Sections 3.0 and 4,0. 
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**«*********p ARAMETRIC flow modulating fan 


*M*0*D*F*A*N* 


AXIAL OR CENT COMPRESSOR? 1-CENT;0 OR CR-AXIAL 
-0 


DESCRIPTION OF INPUT VARIABLES IN MODFAN PROGRAM 
(NAMELIST INPUT) 

ITYPE-1 VIGV 

ITYPE-2 VPF ( NO IGV > 

DESIGN POINT VALUES OF: 


PQPD PRESSURE RATIO! TOTAL -TO-TOTAL ) 

EFFD ADIABATIC EFF IC I ENCY< TOTAL-TO TOTAL) 

W1D INLET CORRECTED FLOW 

PSID ROTOR EXIT MERIDIONAL LOADING! 2GJ*DH/U2**2 ) 

WQA1D INLET CORRECTED FLOW PER UNIT AREA 

XMZ2D MERIDIONAL MACH NO AT ROTOR EXIT 

SMO CONSTANT SPEED STALL MARG I N( PERCENT ) 

STPID IGV METAL ANGLE 

GEOMETRY SPECIFICATIONS 


A3QA2 STATOR INLET TO ROTOR EXIT AREA RATIO 

RRAT INLET HUB TO TIP RADIUS RATIO 

ROQI RATIO OF ROTOR EXIT TO INLET RADII (MERIDIONAL LINE) 

MAP RANGE (ARRAY VALUES MUST BE IN ASCENOING ORDER) 


NR 

AR 

NSPDS 

APCNC 

NSTP1 

ASTP1 

NROT1 

AROT1 


NAMELIST 
ITYPE - 
PQPD - 
W1D 

WQA1D - 
SMD 

A3QA2 - 
ROQI - 
NR 

AR< I ) 

1 

5 

9 


NO OF R VALUES TO BE CALCULATED 
ARRAY OF R VALUESf 1-STALL , 2«MIN-LOSS > 

NO OF SPEED LINES 

ARRAY OF SPEEDS (DECIMAL PERCENT) 

NO OF IGV SETTINGS 

ARRAY OF IGV ANGLESt-10 TO 40 DEG) 

NO OF ROTOR PITCH SETTINGS 

ARRAY OF ROTOR PITCH ANGLES! -5 TO 5 DEG) 


INPUT 

1 . 


1.6550, 

EFFD f 

0.8340, 

223.8000, 

PSID - 

0.8226, 

41 .0700, 

XMZ2D - 

0.4500, 

14.1000, 

STPID « 

0. 

0.9600, 

RRAT * 

0.5000, 

1 .0000, 
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1 .0000, 

1 . 2000 

, 1.4000, 

1 .8000, 

2.0000 

. 2.2000, 

2.6000, 

2 . 6000 

, 3.0000, 


1 . 6000, 
2.4000, 



ORIGINAL PACi* - 
OF POOR QOAU k V 


NSPDS - 6 

APCNC(I), • 

l 0.5000, 0.7000, 0.8000, / 0.9000, 

5 1.0000, 1.1000, 

NSTP1 - 3 

ASTP 1(1 > - 

1 0. 20.0000, 40.0000, 

NROT1 - 0 

AROT1 ( I ) - 

1 0 . 

END NAMELIST INPUT 

ENTER CHANGES TO NAMELIST INPUT 
•SINPUTS 


NAMELIST 
ITYPE - 
POP D - 

INPUT 

1, 

1.6550, 

EFFD 

m 

0.8340 

W1D 

M 

223.8000, 

PSID 

m 

0.8226 

WQA1D 

m 

41 .0700, 

XMZ2D 

m 

0.4500 

SMD 

m 

14. 1000, 

STP1D 

m 

0. 

A3QA2 

m 

0.9600, 

RRAT 

m 

0.5000 

ROQI 

NR 

AR( I ) 

m 

m 

1 .0000, 
11 

m 





1 

1 .0000, 

1.2000, 

1.4000, 

1.6000, 

5 

1.8000, 

2 . 0000 , 

2 . 2000, 

2 . 4000 , 

9 

NSPDS - 

2.6000, 

6 

2.8000, 

3.0000, 


APCNC(I) - 

1 

5 

NSTP1 • 

0.5000, 
1 .0000, 
3 

0.7000, 
1 . 1000, 

0.8000, 

0.9000, 

ASTP 1 < I ) - 

1 

B. 

20.0000, 

40.0000, 



NROT1 - 0 

AROTK I ) - 

1 0 . 

END NAMELIST INPUT 

OUTPUT ON IFC-30,31,32 

PROGRAM COMPLETE & 
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1001 

P-FAN FLOW VS. R. 

SPKFn, AMD 

ANOL 




STP 1 

3 

0 . 0 

20. OOO 

40 . nnu 





SPED 

6 

0 . bOO 

0. 700 

0. 800 

0 . 900 

1 . OOO 

1.100 


R 

1 1 

1 OOO 

1 . 200 

1 . 41 *0 

1 0 . n 

1 . 800 

2 . OOO 

2. 200 

R 

H 

“ “ 2 tOO 

■Kffill!] 

i: . OOO 





FLOW 

1 1 

93. 319/ 

1O0. 9060 

♦ Of. . 071 9 

115.1715 

121 .GO 18 

127. 50 3G 

1 32 . 6690 

FLOW 

1 1 

136. 8715 

1 fi;l . iiflU7 

1 u o007 

1 38 . b 30/ 




FLOW 

1 1 

140 0Mb 

147.4998 

1 r “ i . V ,.>/0 

ILO . v**' 1 

166. 5290 

1 71 . 3347 

1 74 . 6529 

TIThT 

1 1 

Ilf I1TTT1 

T75TTT55 

1 77T TTiTo 

1 /*; y 1 0 




FLOW 

1 1 

164.2152 

1 70. L 422 

1 A* OG 06 

* 3 2:0-3 

186. 8374 

190. 7C5? 

193. 3546 

FLOW 

1 1 

194.0171 

104 . 0171 

194.01 7 1 

19 1.0171 




FL OW 

1 1 

1 87.6732 

192. Pi. 20 

19/ 5544 

i O 1.71/4 

. 32 1 0 

203. 8336 

210. 3036 

FLOW 

1 1 

210. 7646 

210 7010 

2 1 0 /CA~, 

210. 71 : 16 




FLOW 

1 1 

209. 7717 

213. 21 7 0 

'* 6 . 35- •*> 

>13. 1 763 

221 . 6640 

273 . 8120 

225. 3357 

FLOW 

1 1 

225. 9450 

225 94 oO 

-23 V| in 

225. 9440 




FLOW 

1 1 

227. 1011 

229. 7C i7 

231 8 72 3 

7 .3. 7028 

235 5316 

237. 0657 

238. 2C61 

FLOW 

1 1 

238.9715 

P3&. is?4 

.' Jf.. 1 3 24 

2J9 1321 




SPLD 

6 

0 500 

0. 700 

0 . ao- » 

0 . 300 

1 . OOO 

1.100 


R 

1 1 

1 . 000 

1 . 200 

1 .400 

1 . 600 

1 . 800 

2.000 

2.200 

R 

1 1 

2.400 

2. 600 

2 8tM 

3. 0 30 




FLOW 

It ' 

79. 0908 

85. 2.2^5 

?rw 

97. 8971 

1 63 53-1 1 

ion . 7000 

11 503d 

FLOW 

1 1 

117. SwOO 

120.0125 

li. i . 672t 

121 . 873b 




FLOW 

1 1 

117.5 102 

124 . b625 

Cl. 2406 

137 2343 

142.7091 

147.4722 

151 . 2191 

FLOW 

1 1 

103 5704 

154 . 4Mi2 

1*1 476 1 

1 51 4761 




FLOW 

1 1 

137.4297 

1-14. P.TisT 

15*L45/4 

1 5C 05 1 4 

160. 9605 

165 1433 

166. 1275 

FLOW 

1 1 

169.4122 

169.101 1 

)ru, 461 1 

10b. 4Gi 1 




FLOW 

1 1 

157. T r 4i2 

169. 5344 

1C<i 0203 

1/3 80 30 

178.0440 

1/1.5134 

1 S3. 7331 

FLOW 

1 1 

104.2108 

1 34 .2108 

1 r 1.21 06 

1 3 1 21r*3 




FLOW 

1 1 

176. 8737 

181 0057 

180 4038 

1.0. 3422 

193 6995 

1 OG. 45G5 

198. 1221 

FLOW 

1 1 

1 90. 3005 

190. 3605 

193. 3605 

1 9ii . 3b05 




FLOW 

1 1 

195.7194 

199. 3363 

202 bo 7 7 

2i»5 4047 

207.8354 

209.8586 

211.1310 

FLOW 

1 1 

211. 3. 33 

211. 3/33 

21 1 . 37 3*J 

2 l 1 .3733 




SPED 

6 

0 500 

0 . 700 

0 . mo 

0. 900 

1 .OOO 

1.100 


R 

1 1 

1 . uOO 

1 200 

1 . 400 

1 . TOO 

1 . 800 

2. OOO 

2. 200 

R 

1 1 

2.400 

2.000 

2.800 

3.000 




FLOW 

1 1 

67. 3477 

7? 8335 

v Vo. 1246 

83. 1644 

37. 9807 

O'*. 4795 

96.6033 

FI OW 

11 

10O. 2 /09 

103.4467 

1 05 4 11 1 

105.4M1 




FLOW 

1 1 

98 . 8977 

10 b. 09 14 

IIO. 4063 

116.2* MO 

121-1199 

1 :• 5.4542 

1 29 . 0222 

FLOW 

1 1 

131 . 5602 

133.0700 

1 0 98 96 

1 33 4 8*30 




FLOW 

1 1 

1 15. 15,10 

1 2 1 . ,-962 

! i..<‘bO*'l 

|3 ‘ n#»i,2 

136.61 r i 

110.5457 

143.5455 

Ft OW 

1 1 

145 2C62 

145 >835 

1 15 C. 3f 

1 lb. car > 




FLOW 

1 1 

131 .2 *71 

137. 0/71 

1. 12 3G1> 

117. 0791 

151.1319 

Ki 1 . 5567 

1 56. 9?38 

FLOW 

1 1 

157. 7993 

157 763 3 

1 6 / . 7093 

157. 7 . 33 





1 1 

■EHBm 1 ■ 

1»2. 3177 

1 A *.9921 

101.0713 

161. 5343 

167. 3P03 

169.1 323 

FLOW 

1 1 

169.4217 

1 60 4217 

10 » 4217 

160.421 7 




FLOW 

1 1 

162. 3147 

166. 70^5 

I 70 . * **57 

1 73. 6^30 

176. 6320 

1/3.9214 

180. 17 17 

FLOW 

1 1 

180. 2039 

100. 2639 

1 >50 2b 3b 

1 60. 2039 




mm 
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ZZ'Z:. ''4 

OF POOR C 


1002 P-FAN EFF VS. R, SPi r O, AND ANGL 


STP1 

SPED 

R 

3 
6 
1 1 

0. 0 
0 . 500 
1 . ooo 

20 . OOO 
0. 700 
1 . . 00 

4 0 OOO 
0. 800 
1 . ioo 

0 . 900 
1 . oOO 

1 . OOO 
1 . 800 

1.100 
2. OOO 

2. 200 

R 

n — 

2 4 AO 

— 2 pnti 

2 . OOO 

3 . OOO - 




EFF 

1 1 

0. 7 7u<1 

0.8M1 

fi . 6? 1 

O. 8255 

0. 7083 

0.7012 

0 5379 

EFF 

1 1 

0 24-19 

o . ooor 

U . 0006 

0 0006 




EFF 

1 1 

o. n. *rio 

0. 'Ill 

n 8411 

0 0*1 1 

0. 8302 

0.7350 

0. 7100 

FFF 

1 1 

0. 500 2 

0 4273 

0 . 20 70 

0. 00**3 




EFF 

1 1 

o eou 5 

0 0515 

O. Of; 92 

O. 858.0 

0. C460 

0. 8205 

0. 7780 

EFF 

1 1 

0.7178 

0.6502 

f i . VIGO 

0. 430 3 




EFF 

1 1 

0 . 84 1 8 

0 05 15 

n C i 70 

0 8 >7 7 

0.8529 

0.0416 

0. 0231 

EFF 

1 1 

0.7976 

0.7723 

u . 74 1 0 

0. 7025 




EFF 

1 1 

0.8300 

0. 8350 

O. SCSI 

0 . 8390 

0. 837/ 

0. 8340 

0. 8277 

EFF 

1 1 

0.8186 

0.8101 

0.6011 

0.7913 




EFF 

1 1 

O. 7960 

0. 75 86 

0. COH4 

0. 81*1 3 

0.8013 

0. 0003 

0. 7904 

EFF 

1 1 

0. 7955 

0.7915 

0. 765? 

0. 782 8 




SPED 

6 

0. 500 

0. 700 

0 . 1.60 

0 600 

1 . ooo 

1.100 


R 

1 1 

1 ooo 

1 .200 

1 .400 

1 . 600 

1 . ooo 

2.000 

2.200 

R 

1 1 

2.4 GO 

2 600 

2 800 

3 . OOU 




EFF 

1 1 

0 7633 

0.8122 

0 8502 

0. 0749 

0. 8735 

0. 6102 

0 7535 

EFF 

1 1 

0 574 1 

0.2127 

0.0149 

0.0149 




EFF 

1 1 

0 . 7997 

0. 8356 

0. 8694 

0. 8002 

0. 882? 

0. 8034 

0.8157 

EFF 

1 1 

0. 7205 

0 . 571 1 

O 3906 





EFF 

11 

0 6106 

0. 3405 

U . 6637 

0.8702 

0.8814 

0.8694 

0. 8380 

EFF 

1 1 

0. 7809 

0. 71 77 

0 6332 

0 . 5060 




EFF 

1 1 

0.8140 

0. 8363 

0. 0565 

0. 8082 

0.8719 

0.8655 

0 8470 

EFF 

1 1 

0.8190 

0 7904 

O. 7500 

O 6955 




EFF 

1 1 

0 8105 

0.8275 

0.0408 

0.0197 

0.8535 

0.6511 

0.8419 

EFF 

1 1 

0. 8299 

0.8160 

O . 60 1 7 

0 7802 




EFF 

1 1 

0. 7963 

0 . 806 7 

0.6116 

0. 82C5 

0.8234 

0.8231 

0. 8196 

EFF 

1 1 

0.0148 

0.8116 

EkSBB 

0 800 7 




SPED 

6 

0 500 

0. 700 

0 600 

0.900 

1 . OOO 

1.100 


R 

1 1 

1 . OOO 

1 200 

1 .400 

1 . 600 

1 . 800 

2.000 

2. 200 

R 

1 1 

2.400 

2 OOO 

2.600 

a . non 




EFF 

1 1 

0. 7217 

0. 7736 

0 8174 

0.651 1 

0.0681 

0.6595 

0. 0098 

EFF 

1 1 

0 6874 

0. 1155 

0 0126 

0.01 70 




EFF 

1 1 

O. 7491 

0 7901 

O. 8240 

0.8510 

0 . 8G4 7 

0. 0605 

0. 0305 

EFF 

1 1 

0. 7626 

0.0316 

0.4393 

0 . 2 m 




EFF 

1 1 

O 7554 

0. 7910 

0 8 .319 

n. o V.i3 

O. 0533 

0. 8528 

0.6308 

EFF 

1 1 

O 7826 

0.7.34 

0.6522 

0.5511 




EFF 

1 1 

0 . 7555 

0. 7057 

0 8110 

0 82 39 

0. 8102 

0. 0300 

0. 8237 

EFF 

1 1 

0 . 7924 

0. /709 

0 73/3 

0 0830 




EFF 

1 1 

0.7-19? 

0.7796 _ 

O. 7599 

0. 8090 

0.8173 

0.8170 

0.8071 

EFF 

1 1 

0. 7948 

0. 7099 

O. 7066 

0. 7409 




EFF 

1 1 

0. 7350 

0. 7534 

0 . V 065 

0. 7793 

0. 7056 

0.7846 

0. 7796 

EFF 

1 1 

0. 7754 

0. 7707 

0. 7626 

0 . 7509 





£01 
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1003 P TAN PR VS. R, SPEED, AMO ANC L 


STPI 

SPED 

R 

3 
6 
1 1 

0. O 
0. 500 
1 . ono 

20. OOO 
0 700 
1 . 200 

•3. noo 
0 . GOO 
1 . 400 

0.900 
1 . 600 

1 .OOO 
1 . 300 

1.100 
7 . OOO 

2 200 

R 

1 1 

2 400 

2 600 

2 . 80* * 

3 . OOO 




PR 

1 1 

1 . 1507 

1 . 1486 

1 . 1 304 

1 . 1 70 1 

1 . 0990 

1 . 0753 

1 . 0462 

PR 

1 1 

1.01 75 

1 . oooo 

1 . OOOO 

1 . OO. *o 




PR 

1 1 

1 . 323/ 

1 . 3u 70 

1 . 26. > 1 

i . ;\>u/ 

1 . 2275 

1 . 19P0 

1 . 1 526 

PR 

1 1 

1 . 1097 

1 06/3 

1 . 0C*G3 

1 . 0003 




PR 

1 1 

1 . 4007 

1 .1197 

1 . 3951 

1 . 3 .6 1 

1 . 3329 

1 . 2958 

1 . 2551 

PR 

1 1 

1.2119 

1.1714 

1 . 1 304 

1 . 06 3 




PR 

1 1 

1 . 5694 

1 . 56- 5 

1 . ? 11 1 

1 . 51 C 2 

1 .4650 

1 . 4507 

1.4136 

PR 

1 1 

1 . 3745 

1 . 3300 

1 . 300 3 

1 . 2625 




PR 

1 1 

1 . 7698 

1.7510 

1 . 7301 

1 . 70 , 0 

1 . 6820 

1 . C550 

1 . 6262 

PR 

1 1 

1 . 5958 

1 . 5667 

1 . 53 CJ 

1 . 5007 




PR 

1 1 

1 . 9818 

1 . 9664 

1 , 9 107 

1.9316 

1.9127 

1 . 8924 

1 . 8709 

PR 

1 1 

1 . 8485 

1 . 8251 

1 . 804 1 

1 . 7832 




SPED 

6 

0.500 

O. 700 

0. GOO 

0 . 900 

1 . OOO 

1.100 


R 

1 1 

1 .000 

1 . 200 

1 . 400 

1 . 600 

1 . 800 

2.000 

2.200 

R 

1 1 

2.460 

2 GOO 

2.600 

3 . OOO 




PR 

1 1 

1 . 1487 

1 . 1 430 

1 . 1336 

1 . 1 207 

1 . 1045 

1 . 0651 

1 . 0627 

PR 

1 1 

1 . 0376 

1.0104 

l . 0006 

1 . 0006 




PR 

1 1 

1 . 3000 

1 . .864 

1 . 2080 

1 . 2451 

1.2178 

1 . 1866 

1 . 1520 

PR 

1 1 

1.1M6 

1 . 0755 

1 . O il 6 

1 , 0090 




PR 

1 1 

1 . 3906 

1.3017 

1 . 3588 

1.3313 

1 . 2994 

1 . 2631 

1 . 2245 

PR 

1 1 

1 . 1 836 

1 . 1467 

1.1105 

1 . 0736 




PR 

1 1 

1.5162 

1 .4944 

1 .4681 

1 .43 74 

1 ,4026 

1 . 36 42 

1 . 3234 

PR 

1 1 

1 . 2833 

1 . 2462 

1 . 2084 

1 . 1 /no 




PR 

1 1 

1 .65 10 

1 . 6? 38 

1.0017 

1 . 5700 

1 . 5350 

1 .4908 

1 .4570 

PR 

1 1 

1.4192 

1 . 3824 

1 . 3470 

1.3101 




PR 

1 1 

1 . 0091 

1 . 7e 19 

1 . 7500 

1 . 72 r 3 

1 . 6969 

1 . 6631 

1 . 6284 

PR 

1 1 

1 . 5955 

1 . 5650 

1 . 5342 

1 . 5029 




SPED 

6 

0.500 

0.700 

0. 800 

0 . 900 

1 . OOO 

1.100 


R 

1 1 

1 . oon 

1 200 

1 400 

1 . 600 

1 . 800 

2.000 

2 200 

R 

1 1 

2.400 

2 . 600 

2 . 000 

3.000 




PR 

1 1 

1 . 1 360 

1.1310 

1 . 123U 

1.1121 

1 . 0682 

1 . 0316 

1 . 0626 

PR 

1 1 

1.0414 

1.0182 

1 . 0004 

1 . 0004 




PR 

1 1 

1 . 2 < 37 

1 . 2577 

1 . 241* 

1 . 2203 

1 . 1 965 

1 . 1687 

1.1301 

PR 

1 1 

1 . 1059 

1 .0721 

1 . 0299 

1.0150 




PK 

11 

1.3528 

1 . 3363 

1.3152 

1 . 2697 

1 . 2503 

1 . 2273 

1 . 1 920 

PR 

1 1 

1 . 1556 

1 . 1220 

1 . 0938 

1 . 0047 




PR 

1 1 

1 . 4470 

1 . 4250 

1 .4000 

1 . 36M9 

1 . 3339 

1 . 2966 

1 . 2596 

PR 

1 1 

1.2214 

1 . 1 90b 

1 . 1 509 

1 . 1 261 




PR 

1 1 


1 .W10 

1 _1_951 

1 .4613 

1 . 4241 

1 . 3839 

1 . 3431 

PR 

1 1 

1 . 3074 

1 . 2745 

1 . 2409 

1 . 2066 




PR 

1 1 

1 . 6604 

1.6315 

1 . 5u»91 

1 . b» 33 

1 . 5248 

1 .4838 

1 . 4433 

PR 

1 1 

1 . 4096 

1 . 3767 

1 . 34 32 

1 . 3091 
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Ouidi ku i - -' - ' -• 

OF POOR QUALhV 


PARAMETRIC FLOW MODULATING FAN 


*M*0*D*F*A*N* 


AXIAL OR CENT COMPRESSOR? 1-CENT;0 OR CR-AXIAL 

-I 


DESCRIPTION OF INPUT VARIABLES IN MODFAN PROGRAM 
(NAMELIST INPUT) 


ITYPE-1 V IGV 

ITYPE »2 VPF ( NO IGV ) 


DESIGN POINT VALUES OF: 


PQPD 

EFFD 

W1D 

PSID 

VQA1D 

XM22D 

SMO 

STP1D 


PRESSURE RATIO< TOTAL-TO-TOTAL > 

ADIABATIC EFFICIENCY( TOTAL -TO TOTAL) 

INLET CORRECTED FLOW 

ROTOR EXIT MERIDIONAL LOADING( 2GU*DH/U2**2 ) 
INLET CORRECTED FLOW PER UNIT AREA 
MERIDIONAL MACH NO AT ROTOR EXIT 
CONSTANT SPEED STALL MARG I N < PE RCENT ) 

IGV METAL ANGLE 


GEOMETRY SPECIFICATIONS 


A3QA2 

RRAT 

ROQI 


STATOR INLET TO ROTOR EXIT AREA RATIO 
INLET HUB TO TIP RADIUS RATIO 

RATIO OF ROTOR EXIT TO INLET RADII (MERIDIONAL LINE) 


MAP RANGE (ARRAY VALUES MUST BE IN ASCENDING ORDER) 


NR 

AR 

NSPOS 

APCNC 

NSTP1 

ASTP1 

NROT1 

AROT1 


NO OF R VALUES TO BE CALCULATED 
ARRAY OF R VALUE S< 1 -STAL L , 2-MI N-LOSS ) 

NO OF SPEED LINES 

ARRAY OF SPEEDS (DECIMAL PERCENT) 

NO OF IGV SETTINGS 

ARRAY OF IGV ANGLES(-10 TO 40 DEG) 

NO OF ROTOR PITCH SETTINGS 

ARRAY OF ROTOR PITCH ANGLES( -5 TO 5 DEG) 


NAMELIST INPUT 




ITYPE - 

1 . 




PQPD - 

4,3000, 

EFFD - 

0.7600, 


W1D 

5.5000, 

PSID - 

1 .8160, 


WQA1D - 

35.3100, 

XMZ2D - 

0.2950, 


SMD 

10.0000, 

STP1D « 

0. 


A3QA2 « 

1 ,0500, 

RRAT - 

0.2750, 


ROQI - 

2.0300, 




NR 

1 1 




AR< l ) 

• 




1 

1 .0000, 

1 . 2000 

, 1.4000, 

1 . 6000 

5 

1 .8000, 

2.0000 

, 2.2000, 

2 . 4000 

9 

2.6000, 

2 . 8000 

, 3.0000, 




NSPOS - 6 

APCNC(I) - 

1 0.5000, 0.7000, 0.8000, 0.9000, 

S 1.0000, 1.1000, 

NSTP1 - 1 

ASTPl(I) - 

1 0 . 

NROT1 - 0 

AROT1 < I ) - 

1 0 . 

END NAMELIST INPUT 

ENTER CHANGES TO NAMELIST INPUT 
-SINPUTS 

NAMELIST INPUT 

ITYPE - 1, 


PQP D 

m 

4 . 3000 , 

EFFD 

■ 

0.7600, 

W1D 

s 

5.5000, 

PSID 

a 

1.8160, 

WQA10 

a 

35.3100, 

XMZ2D 

■ 

0.2950, 

SMD 

m 

10.0000, 

STP10 

a 

0. 

A3QA2 

■ 

1 .0500, 

RRAT 

a 

0.2750, 

ROQI 

a 

2.0300, 





NR 11 

AR( I ) 

1 1.0000, 1.2000, 1.4000, 1.6000, 

5 1.8000, 2.0000, 2.2000, 2.4000, 

9 2.6000, 2.8000, 3.0000, 

NSPDS - 6 

APCNCII) - 

1 0.5000, 0.7000, 0.8000, 0.9000, 

5 1.0000, 1.1000, 

NSTP1 - 1 

ASTPH I ) * 

1 0 . 

NROT1 - 0 

AROT1 ( I ) - 

1 0 . 

ENO NAMELIST INPUT 

OUTPUT ON IFC-30,31,32 

PROGRAM COMPLETE & 
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c-r: ; ;' " " r 

OF r'Gi n v/ 


1001 


P-FAN FLOW 

VS. R, 

SPEED i AND 

ANGL 




STP1 

1 

0 0 







SPED 

6 

0. 500 

0. 700 

0. noo 

0 . 900 

1 . 000 

1.100 


R 

1 1 

1 . GOO 

1 . 200 

1 . 400 

1 . GOO 

1 . 0OO 

2.000 

2 200 

R 

1 1 

2 400 

2.600 

2 . 000 

3 . 000 




FLOW 

1 1 

2. 1 404 

2 2372 

2. 2270 

2. 4156 

2 . 5029 

2. 5860 

2. 6730 

FLOW 

1 1 

2. 75ii3 

2. 0355 

P . 9 1 35 

2 939 1 




FLOW 

1 1 

3. 20/0 

3. 3132 

3.41/5 

3.5198 

3. 0200 

3.7100 

3. 6127 

FLOW 

1 1 

3 9032 

3. 0092 

4 . 0/06 

4. 14G9 




FLOW 

1 1 

3 . 79 1 2 

3. 0982 

4 . 0033 

4. 1062 

4 . 2069 

4 . 3051 

4.3997 

FLOW 

1 1 

4.4892 

4.5732 

4 6517 

4 . 7213 




FLOW 

1 1 

4.4077 

4.5106 

4 . G 1 1 G 

4. 7105 

4.8073 

4 . 901 7 

4 . 9924 

FLOW 

1 1 

5.0777 

5. 1 572 

5. 2308 

5.2902 




FLOW 

1 1 

5. 0481 

5. 1420 

5. 23 42 

5. 3246 

5.4132 

5.4999 

5. 5831 

FLOW 

1 1 

5.661 1 

5.7337 

5. 0008 

5 . 3620 




FLOW 

1 1 

5.6978 

5. /771 

5 0552 

5 . 9320 

6. 0075 

6.0815 

6 . 1 527 

FLOW 

1 1 

6.2197 

6. 2824 

G. 3405 

6. 3940 





EOT 
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1002 P-FAN EFF VS. R, SPEED, AND ANGL 

STP 1 1 0.0 


SPED 

6 

0.500 

0.700 

0.800 

0 900 

1 . 000 

1.100 


R 

1 1 

1 . 000 

1 . 200 

i . ino 

1 . 600 

1 . 800 

2.000 

2 200 

R 

1 1 

2 . *100 

2. 600 

2 . 800 

3 OOO 




EFF 

1 1 

0.8195 

0.8210 

0. 621 7 

0.8313 

0.0199 

0.8175 

0.8142 

EFF 

1 1 

0. 8098 

0. 8045 

0. 7082 

0 . 7909 




EFF 

1 1 

0. 0259 

0. f.2 72 

0. 8273 

0 . 8278 

0. 8272 

o. a^>r>o 

0. 8242 

EFF 

1 1 

0.0217 

0.8187 

0. 8150 

0. 8107 




EFF 

1 1 

0.8136 

0.8147 

0.815? 

0.81 64 

0.8151 

0 8143 

0.8130 

EFF 

1 1 

0.8113 

0. 8002 

0. 0006 

0 8035 




EFF 

1 1 

0.7918 

0. 7026 

0 . 7 932 

0. / 933 

0. 7932 

0. 7920 

0. 7920 

EFF 

1 1 

0. 7909 

0. 7C95 

0. 7078 

0. 7850 




EFF 

1 1 

0 7590 

0 . 7596 

0. 7600 

0. 7t,*)2 

0. 7602 

0. 7600 

0. 7596 

EFF 

1 1 

0. 7589 

0. 75 01 

0 . / 5 / 1 

0.7559 




EFF 

1 1 

0. 7103 

0.7107 

0.7110 

0.7111 

0.7112 

0.7111 

0. 7109 

EFF 

1 1 

0.7106 

0.7101 

0 . 7096 

0. 7089 





EOT 
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OF POOR QUALITY 


1003 P-FAN PR VS. R, SPEED, AND ANGL 

STP1 1 0.0 


SPED 6 

0.500 

0. 700 

0. 800 

0.900 

1 .OOO 

1.100 


R 1 1 

1 . OOO 

1 . 200 

1 . 400 

1 . 600 

1 . 500 

2 OOO 

2.200 

R 1 1 

2 . 400 

2 600 

2. 800 

3 . OOO 




PR 1 1 

1 . 5706 

1 . 5749 

1 . 5723 

1 . 81.90 

1 . 5650 

1 . 5601 

1 5544 

PR 1 1 

1 . 5451 

1.5410 

1 . 5333 

1 5218 




PR 1 1 

2. 3420 

2. 3367 

2. 33 >0 

2. 3731 

2. 3708 

2.3123 

2. 3025 

PR 1 1 

2.2915 

2. 2734 

2.2601 

2.2517 




PR 1 1 

2. 8901 

2. 86 ‘0 

2. 6864 

2. 6793 

2. 0708 

2. 8607 

2. 8193 

PR 1 1 

2. 8305 

2 823 1 

2. 8070 

2 . 7903 




PR 1 1 

3. 5714 

3. 50 »7 

3 . 5005 

3. 85>’C 

3.5435 

3.5327 

3 . 5704 

PR 1 1 

3. 5066 

3.4918 

3 4754 

3.4576 




PR 1 1 

4. 3386 

4 . 3337 

4 . 3274 

4.3197 

4.3105 

4 . 3000 

4. 2880 

PR 1 1 

4.2748 

4 . 2602 

4 . 2144 

4 . 2274 




PR 1 1 

5.0364 

5.0616 

5 . 0 758 

5 . 0600 

5 . 0606 

5. 051 3 

5 04 10 

PR 1 1 

5.0296 

5.0171 

5.0036 

4 . 9891 





EOT 



7.0 ANALYTICAL BACKGROUND 


7.1 Method of Map Generation 

The fan geometry is defined implicitly by the design point input, and the user 
selected option such as VIGV settings etc. The rotor design parameters are first 
separated from those of the IGV and stator. The rotor is then analyzed separately and 
the IGV and stator losses added on after the rotor calculation has been completed. 

The min-loss line which forms the backbone of the map is assumed to pass through 
the design point. This sets the optimum incidence angle on the rotor. The flow co- 
efficient (Cz^AJi) is assumed to remain constant with speed along the min-loss line 

2 

for any given value of STP1 or R0T1. The work coefficient (2g J DH/U2 ) is obtained 
from the flow coefficient, rotor geometry, and rotor continuity. The rotor loss at 
min-loss is obtained at the design point. This loss is assumed to be a function of 
rotor inlet relative Mach number only since at design point the incidence is assumed 
optimum. A curve of loss vs. rotor inlet relative Mach number is used with a scalar 
to force the loss through the design point value. For off-design the incidence loss is 
assumed to vary as (l-cos^i) where n is determined from the NASA TASK II data of Ref. 4. 

In order to indicate how different IGV positions are calculated, consider the 
min-loss point on the 100% speed line. Assume an IGV setting of STP1= +20° is required. 
The rotor exit relative flow angle (BETA2) is assumed to remain constant. The flow 
coefficient and work coefficient are then calculated from the known value of ALPRA1 and 
BETA1. A slight shift in the min-loss value of BETA1 and IGV position has been observed 
in the data and is applied as a correction to BETA1 before the calculation is made. 

The rotor continuity iteration to determine the axial valocity ratio across the 
rotor is the key iteration in the program. The results of this iteration together with 
the known flow angles sets the work coefficient at min-loss. 

To determine the off-backbone characteristics, the angles ALPHA1 and BETA2 are 
assumed equal to their min-loss values at the selected speed. A value of the work 
coefficient is then chosen and the rotor continuity iteration carried out. The stator 
and IGV losses are then added to obtain the stage performance. 
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7.2 


•-■r v ' 

Discussion of Variable Geometry Options for Axial Flow Fans 

7.2.1 Variable Inlet Guide Vane (VIGV) Option 

The VIGV option can be used to generate a set of off-design performance maps 
for a user selected design point and set of IGV angles (i.e., STP1) . The design 
point is assumed to be at a zero value of STP1. The VIGV are assumed to be of the 
type tested in the NASA TASK II program (Ref. 4), i.e., of a flap-type as sketched 
in Figure 4. Since the IGV leading edge does not move, no IGV incidence loss is in- 
cluded in the stage loss calculation. The inlet flow angle (ALPHA1) is somewhat less 
than the IGV angle (STP1) . As currently set the ratio of ALPHAl to STP1 is about 
0.775, as determined from the data of Ref. 4. 

The nature of the flow modulation produced by the VIGV’s is illustrated in 
Figure 4. Let the nominal IGV setting represent the design point. Then at the same 
speed and rotor incidence angle, the closed position results in a smaller value of 
inlet axial velocity and, therefore, a smaller flow. A similar line of reasoning 
leads to the conclusions that a negative IGV setting results in a flow increase. 

The above conclusions can be made somewhat clearer and more quantative by 
referring to the stage characteristic sketch shown in Figure 5. The fan stage 
characteristic can be expressed in either of the following two ways. 


t = 2-Z <t>, + ^ , fc c (Sj 0) 

(4 - ^ ( fa 0> t fc* 0>J ( 2 ) 

where $4= LH/v£/i< t T * <p, * 

d, * absolute Inlet angle* 

■ relative inlet angle 

relative exit angle 
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Figure 4. Vector Diagram for Variable Inlet Guide Vanes (VIGV) 
at Constant Rotor Incidence and Deviation Angle. 
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Figure 5. Stage Characteristics for VIGV Axial-Flow Fan. 
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If the axial velocity ratio (C z2 /C z i) and the exit flow angle are assumed to be con- 
stant, then equation one is a straight line passing through two (2) when the flow 
coefficient equals zero. A series of these lines for different values of STP1 
(ALPHA1 = 0.775 x STP1) are shown in Figure 5. Also shown in Figure 5 is a line of 
constant incidence angle passing through the design point (as given by equation 2) . ' 
The locus of the min— loss point on the 100% speed line, as generated by the programs, 
has also been indicated in the figure. It can be seen from Figure 5 that the min-loss 
locus closely follows a line of constant rotor incidence angle. The changing flow 
coefficient explains the flow modulation observed. Parametric maps generated for 
the NASA TASK II VIGV Comparison (Ref. 4) are shown in Figures 6 thru 9. One negative 
angle map was generated from the program in order to show the basic trends. Test data 
is reported in (Ref. 4) includes three positive IGV positions (0°, 20°, and 40°) and 
no negative IGV settings. 

The test data of Reference 4 was compared with the program results along the 
min-loss locus. These comparison plots are shown in Figure 10, 11 and 12. As can be 
seen from the plots, as the IGV are closed, both the flow and the pressure rise de- 
crease. In general, the agreement between the map and the test data is quite good. 

All of the trends in the data are correctly predicted by the program. 

7.2.2 Variable Pitch Rotor (VPF) Option 

The VPF option can be used to generate a set of off-design performance maps 
for a user selected design point and set of rotor pitch angles (i.e., R0T1) . As with 
the VIGV option, the design point is assumed to be at a zero value of R0T1. IGV are 
assumed to be absent in this option. 

The nature of the flow modulation produced by the VPF is illustrated in Figure 
13. Let the nominal pitch setting represent the design point. Then at the same speed 
and zero incidence angle the positive pitch (i.e., rotor closing) has a smaller value 
of axial velocity than the nominal pitch. This results in less flow as the rotor is 
closed. 

As with VIGV, the same conclusions follow from the stage characteristic. This 
is shown in Figure 14. Note thate the min-loss locus generated by the program closely 
follows zero rotor incidence angle. The relative flow angle corresponding to zero 
incidence is changed by varying the pitch of the rotor. Note that the min-loss locus 
is nearly horizontal for the VPF while for the VIGV's the locus sloped upward from 
the origin. 
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Figure 6. Parametric IVGV Axial Flow Fan (STP1 = -10.0°). 
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Figure 11. Pressure Rise Variation Along Min-Loss Locus (Single 
Stage VIGV Fan) . 
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Figure 13. Vector Diagram for Variable Pitch Rotor (VPF) . 
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7.3 


Discussion of Centrifugal Compressor Option 


7.3.1 Fixed Geometry Centrifugal Compressors 

The centrifugal compressor option can be used to generate a set of off-design 
performance maps for a user selected design point. The default settings assume a 
geometry consisting of an impellor followed by a radial diffuser. The inlet swirl 
angle is assumed to be zero and no provision is made for inlet guide vane loss. 


The test data of Reference 5 was compared with the program results for the 
centrifugal compressor default case along the min-loss locus. These comparison plots 
are shown in Figure 15, 16, and 17. In general, the agreement between the map and the 
test data is fairly good. All of the trends in the data are correctly predicted by the 
program. 


7.3.2 Variable Inlet Guide Vanes (VIGV) 

Centrifugal compressors frequently have a substantial degree of inlet swirl. 
This swirl may result from the presence of an upstream axial compressor, or from the 
presence of inlet guide vanes. In order to account for inlet swirl at the design point, 
the program has provision for user selection of the design point inlet guide vane 
position (STP1D) . The turning effectiveness of the IGV is assumed to 0.775 (i.e., 

ANGA1 = 0.775 x STP1) . No loss is charged to the IGV. 

A variable inlet guide vane option (VIGV) has been provided for the centri- 
fugal compressor just as for the axial flow fan. The most significant difference 
between the axial and centrifugal VIGV options is that no loss is charged to the inlet 
guide vanes in the latter case. 


The nature of the flow modulation produced by the VIGV's for the centrifugal 
machine is somewhat different from that for an axial machine. This can be made clear 
by comparing the stage characteristics for the VIGV centrifugal compressor with those 
of the axial flow fan which was shown in Figure 5. 

The stage characteristic for a centrifugal compressor can be expressed in 
either of the following two ways; 


{ I P, Yf fix ><L 

<£/ Z £/?, 




r A 
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Figure 16. Pressure Rise Variation Along Min-Loss Locus 
(Centrifugal Compressor) . 











Equations three and four reduce to equations one and two if the ratio of the inlet 
to exit radii is set equal to unity. 

If the axial velocity ratio (C z 2 /C z i) and the exit flow angle are assumed to 
be constant, then equation three is a straight line passing through two (2) when the 
flow coefficient equals zero. A series of these lines for different values of STP1 
(ALPHA1 = 0.775 x STP1) are shown in Figure 18. Also shown in Figure 18 is a line of 
constant incidence angle passing through the design point (as given by equation 4). 

The locus of the min-loss point on the 100% speed line has also been indicated in the 
figure. It can be seen from Figure 18that the min-loss locus closely follows a line 
of constant rotor incidence angle. The flow modulation results from the change in 
flow coefficient along this locus. 

Note that the value of the min-loss work coefficient increases as the flow coefficient 
is reduced. This behavior contrast with that of the VIGV axial flow fan shown in 
Figure 5 , where a reduction in the value of the min-loss flow coefficient resulted 

in a reduction in the work coefficient. The difference is due to the presence of the 
inlet-to-exit radius ratio in equations three and four. 

Performance maps were generated for the VIGV centrifugal compressor reported 
in Reference 6 . In this reference test data was reported for three IGV positions 
13 , 0 , 40 ). The data as reported gives only the impeller performance. In order to 
generate a map of the impeller performance the program source was edited to eliminate 
the diffusser loss. The test data of Reference 6 was compared with the program results 
along the min-loss locus. These comparison plots are shown in Figures 19, 20, and 21. 

As can be seen from the plots as the IGV are closed, the flow is reduced and the impellei 
pressure rise increases. The agreement between the map and the test data is fairly good. 
According to Reference 6 , an additional loss of about 2.5% in impeller efficiency at 
the STP1 = 40° position occurred as a result of inlet flow distortion. 
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Figure 18. Stage Characteristic for VIGV 
(Centrifugal Compressor) . 
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Figure 20. Pressure Rise Variation Along Min-Loss Locus 
(VIGV Centrifugal Compressor-Impeller Only) . 
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